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Abstract. We present a review of recent work on the statistical mechanics of 
non equihbrium processes based on the analysis of large deviations properties 
of microscopic systems. Stochastic lattice gases are non trivial models of such 
phenomena and can be studied rigorously providing a source of challenging 
mathematical problems. In this way, some principles of wide validity have 
been obtained leading to interesting physical consequences. 



1. A Physicist motivation 

In equilibrium statistical mechanics there is a well defined relationship, estab- 
lished by Boltzmann, between the probability of a state and its entropy. This fact 
was exploited by Einstein to study thermodynamic fluctuations. So far it does 
not exist a theory of irreversible processes of the same generality as equilibrium 
statistical mechanics and presumably it cannot exist. While in equilibrium the 
Gibbs distribution provides all the information and no equation of motion has to 
be solved, the dynamics plays the major role in non equilibrium. 

When we are out of equilibrium, for example in a stationary state of a system in 
contact with two reservoirs, even if the system is in a local equilibrium state so that 
it is possible to define the local thermodynamic variables e.g. density or magneti- 
zation, it is not completely clear how to define the thermodynamic potentials like 
the entropy or the free energy. One possibility, adopting the Boltzmann-Einstein 
point of view, is to use fluctuation theory to define their non equilibrium analogs. 
In fact, in this way extensive functionals can be obtained although not necessarily 
simply additive due to the presence of long range correlations which seem to be a 
rather generic feature of non equilibrium systems. 

Let us recall the Boltzmann-Einstein theory of equilibrium thermodynamic fluc- 
tuations. The main principle is that the probability of a fluctuation in a macroscopic 
region of fixed volume V is 

P cx exp{V^A5/fc} (1.1) 

where AS* is the variation of the specific entropy calculated along a reversible trans- 
formation creating the fluctuation and k is the Boltzmann constant. Eq. Hl.l|l was 
derived by Einstein simply by inverting the Boltzmann relationship between en- 
tropy and probability. He considered Hl.l|l as a phenomenological definition of the 
probability of a state. Einstein theory refers to fluctuations for equilibrium states, 
that is for systems isolated or in contact with reservoirs characterized by the same 
chemical potentials. When in contact with reservoirs AS" is the variation of the 
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total entropy (system + reservoirs) which for fluctuations of constant volume and 
temperature is equal to — AjF/T, that is minus the variation of the free energy of 
the system divided by the temperature. 

We consider a stationary non-equilibrium state (SNS), namely, due to external 
fields and/or different chemical potentials at the boundaries, there is a flow of 
physical quantities, such as heat, electric charge, chemical substances, across the 
system. To start with, it is not always clear that a closed macroscopic dynamical 
description is possible. If the system can be described by a hydrodynamic equation, 
a fact which can be rigorously established in stochastic lattice gases, a reasonable 
goal is to find an explicit connection between the thermodynamic potentials and 
the dynamical macroscopic properties like transport coefficients. The study of large 
fluctuations provides such a connection. 

Besides the definition of thermodynamic potentials, in a dynamical setting a 
typical question one may ask is the following: what is the most probable trajectory 
followed by the system in the spontaneous emergence of a fluctuation or in its 
relaxation to an equilibrium or a stationary state? To answer this question one 
flrst derives a generalization of the Boltzmann-Einstein formula from which the 
most probable trajectory can be calculated by solving a variational principle. For 
equilibrium states and small fluctuations an answer to this type of questions was 
given by Onsager and Machlup in 1953 [231 • The Onsager-Machlup theory gives 
the following result under the assumption of time reversibility of the microscopic 
dynamics: the most probable creation and relaxation trajectories of a fluctuation 
are one the time reversal of the other. 

We discuss this issue in the context of stochastic lattice gases in a box of linear 
size N with birth and death process at the boundary modeling the reservoirs. We 
consider the case when there is only one thermodynamic variable, the local density 
denoted by p. Its macroscopic evolution is given by the continuity equation 



where D(p) is the diffusion matrix, x{p) the mobility and E the external field. 
Here J{p) is the macroscopic instantaneous current associated to the density profile 
p. Finally the interaction with the reservoirs appears as boundary conditions to 
be imposed on solutions of (|1.2|) . We shall denote by u the macroscopic space 
coordinate and by p = p{u) the unique stationary solution of H1.2() . i.e. p is the 
typical density profile for the SNS. 

This equation derives from the underlying stochastic dynamics through an ap- 
propriate scaling limit in which the microscopic time and space coordinates are 
rescaled diffusively. The hydrodynamic equation H1.2|l thus represents the law of 
large numbers for the empirical density of the stochastic lattice gas. The conver- 
gence has to be understood in probability with respect to the law of the stochastic 
lattice gas. Finally, the initial condition for 1)1.2(1 depends on the initial distribution 
of particles. Of course many microscopic configurations give rise to the same initial 
condition pa(u). 

Let us denote by the invariant measure of the stochastic lattice gas. The 
free energy T{p), defined as a functional of the density profile p = p{u), gives 
the asymptotic probability of fluctuations of the empirical measure tt^ under the 
invariant measure . More precisely 



dtp = V- D{p)\/p - x{p)E 



V • J{p) 



(1.2) 




(1.3) 
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where d is the dimensionahty of the system, tt « p means closeness in the weak 
topology and ~ denotes logarithmic equivalence as — > oo. In the above formula 
we omitted the dependence on the temperature since it does not play any role in 
our analysis; we also normalized J- so that J-{p) = 0. 

In the same way, the behavior of space time fluctuations can be described as 
follows. Let us denote by P^n the stationary process of the stochastic lattice gas, 
i.e. the initial distribution is given by the invariant measure . The probability 
that the evolution of the random variable tt^ deviates from the solution of the 
hydrodynamic equation and is close to some trajectory pt is exponentially small 
and of the form 

V,.{7r^ ^pt,te[h,t2]) ^e^pi-N^lTiptJ+Iit.MiP)]} (1-4) 
where I{p) is a functional which vanishes if pt is a solution of (|1.2I) and ^{pti) is the 
free energy cost to produce the initial density profile pt^ . Therefore I{p) represents 
the extra cost necessary to follow the trajectory pt in the time interval [^1,^2]- 

To determine the most probable trajectory followed by the system in the spon- 
taneous creation of a fluctuation, we consider the following physical situation. The 
system is macroscopically in the stationary state p at t = —00 but at t = we find 
it in the state p. According to H1.4(l the most probable trajectory is the one that 
minimizes / among all trajectories pt connecting p to p in the time interval [— cxd, 0], 
that is the optimal path for the variational problem 

F(p)-inf/[_oo,o](p) (1.5) 
p 

The functional V{p), called the quasi-potential, measures the probability of the 
fluctuation p. Moreover, the optimal trajectory for (|1.5|1 determines the path fol- 
lowed by the system in the creation of the fluctuation p. As shown in |P 12 E| this 
minimization problem gives the non equilibrium free energy, i.e. ^ = As we 
discuss here, by analyzing this variational problem for SNS, the Onsager-Machlup 
relationship has to be modified in the following way: the spontaneous emergence of 
a macroscopic fluctuation takes place most likely following a trajectory which can 
be characterized in terms of the time reversed process. 

Beside the density, a very important observable is the current flux. This quantity 
gives informations that cannot be recovered from the density because from a density 
trajectory we can determine the current trajectory only up to a divergence free 
vector field. We emphasize that this is due to the loss of information in the passage 
from the microscopic level to the macroscopic one. 

To discuss the current fluctuations in the context of stochastic lattice gases, we 
introduce the empirical current which measures the local net flow of particles. 
As for the empirical density, it is possible to prove a dynamical large deviations 
principle for the empirical current which is informally stated as follow. Given a 
vector field j : [0, T] x A R'', we have 

P^«(u;^«j(i,^^)) ~exp{-iV'^X[o,T](j)} (1.6) 
where P^n is the law of the stochastic lattice gas with initial condition given by 
— {j^^^-i which represents the number of particles in each site, and the rate 
functional is 

%T](j) = \ f^dt{{2-J{p%x[pr\3-Jm (1.7) 
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in which we recall that 

J{p)^-D{p)Vp + x{p)E . 

Moreover, p = p{t, u) is obtained by solving the continuity equation dtp + S7 ■ j = 
with the initial condition p(0) = po associated to rj^ . The rate functional vanishes if 
j = J{p), so that p solves H1.2|l . This is the law of large numbers for the observable 
. Note that equation 1)1. 7|l can be interpreted, in analogy to the classical Ohm's 
law, as the total energy dissipated in the time interval [0, T] by the extra current 
J-Jip)- 

Among the many problems we can discuss within this theory, we study the 
fluctuations of the time average of the empirical current over a large time interval. 
We show that the probability of observing a time-averaged fluctuation J can be 
described by a functional $(J) which we characterize in terms of a variational 
problem for the functional I[o^t] 

HJ)^\im infiX[o,T](j), (1-8) 

where the infimum is carried over all paths j = u) having time average J. We 
finally analyze the variational problem H1.8() for some stochastic lattice gas models 
and show that different scenarios take place. In particular, for the symmetric exclu- 
sion process with periodic boundary condition the optimal trajectory is constant in 
time. On the other hand for the KMP model 22 , also with periodic boundary con- 
ditions, this is not the case: we show that a current path in the form of a traveling 
wave leads to a higher probability. 



2. Boundary driven simple exclusion process 

For an integer > 1, let A^r := {1, . . . , A^ — 1}. The sites of Aat are denoted by 
a;, 2/, and z while the macroscopic space variable (points in the interval [0,1]) by 
u. We introduce the microscopic state space as Sat := {0, 1}^" which is endowed 
with the discrete topology; elements of Sjv, called configurations, are denoted by 
r\. In this way riix) G {0, 1} stands for the number of particles at site x for the 
configuration r\. 

The one dimensional boundary driven simple exclusion process is the Markov 
process on the state space S^r with infinitesimal generator defined as follows. Given 
a,/3 G (0, 1) we let 

/V2 

(L^/)(,y) := _ ^ [/(a-.^+ir;) - /(r;)] 

+ — - ^(1)} + (1 - «)^(i)] [/(^S) - fiv)] 

+ — [/3{1 - viN 1)} + (1 - PHN 1)] [f{a^-'v) - fiv)] 

for every function / : S^r ^ R. In this formula a^'^r] is the configuration obtained 
from T] by exchanging the occupation variables ri{x) and 77(2/): 

{ri{y) \i z = X 
i]{x) if z = y 
r]{z) ifz^x,y 
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and a^T] is the configuration obtained from rj by flipping the configuration at x: 

where 5x,y is the Kronecker delta. The parameters a, (3, which affect the birth and 
death rates at the two boundaries, represent the densities of the reservoirs. Without 
loss of generality, we assume a < (3. Notice finally that Ln has been speeded up 
by N'^; this corresponds to the diffusive scaling. 

The Markov process {rjt : t > 0} associated to the generator ijv is irreducible. It 
has therefore a unique invariant measure, denoted by i^^ p. The process is reversible 
if and only if a = /3, in which case z/^^, is the Bernoulli product measure with density 
a 

.N 



Va,a{ri ■■ vix) = 1} = a 

ioT 1< X < N - I. 



li a P the process is not reversible and the measure a carries long range 



correlations. Since E^n ^[Lpfj]{x)] = 0, it is not difficult to show that p^{x) 
E^N ^[ri(x)] is the solution of the linear equation 



Ajvp^(x)=0, l<a;<iV-l 
p^(0) = a, p^{N)=f3, 

where Ajv stands for the discrete Laplacian. Hence 



(2.1) 



p''ix)^a + -i/3^a) (2.2) 
It is also possible to obtain a closed expression for the correlations 

K^Jv{x);v{y)] = K^Jv{x)riiy)] - K^Jvix)]E^^j7]{y)] 

As shown in ^2 |2S] , for 1 < a: < j/ < — 1 we have 

E.,JvixyMy)]-~^-^^{i-^) (2.3) 

Computing LNrj{x)rj{y)^ we obtain that the correlations solve a discrete differential 
equation. One can then check that 1)2. 3|l is the solution. 

Note that, if we take x,y at distance 0{N) from the boundary, then the co- 
variance between r](x) and r]{y) is of order 0{\/N). Moreover the random vari- 
ables rj{x) and ■q{y) are negatively correlated. This is the same qualitative behav- 
ior as the one in the canonical Gibbs measure given by the uniform measure on 
SAr,fc = {77 e Sat : YflZi v{x) = k}- 

3. Stationary large deviations of the empirical density 

Denote hy the space of positive measures on [0, 1] with total mass bounded 
by 1. We consider A4+ endowed with the weak topology. For a configuration 77 in 
Sjv, let TT^ be the measure obtained by assigning mass to each particle and 
rescaling space by iV^^ 

x=l 
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where 6u stands for the Dirac measure concentrated on u. Denote by {tt^ , H) the 
integral of a continuous function i7 : [0, 1] ^ M with respect to tt^ 



N-l 



Um 

JV- 



x = l 

We use the same notation for the inner product in L2{[0, l\,du). Analogously we 
denote the space integral of a function / by (/) — du f{u). 

The law of large numbers for the empirical density under the stationary state 
u^jf is proven in [1111111111 ■ 

Theorem 3.1. For every continuous function _ff : [0, 1] — > M and every S > 0, 
where 

p{u) = a{l - u) + l3u . (3.1) 
We remark that p is the solution of the elliptic linear equation 

p(0) = a , p(l) - (3 . 

which is the continuous analog of (|2.1|l . Here and further A stands for the Lapla- 
cian. 

Once a law of large numbers has been established, it is natural to consider 
the deviations around the typical value p. From the explicit expression of the 
microscopic correlations H2.3() it is possible to prove a central limit theorem for the 
empirical density under the stationary measure i^^^- We refer to |25j for a more 
detailed discussion and to for the mathematical details. 

Fix a profile 7 : [0, 1] — > [0, 1] different from p and a neighborhood Ve{'y) of radius 
e > around the measure j{u)du in M + . The mathematical formulation of the 
Boltzmann-Einstein formula Hl.l|l consists in determining the exponential rate of 
decay, as j 00, of 

Derrida, Lebowitz and Speer jl2l I13j derived, by explicit computations, the large 
deviations principle for the empirical density under the stationary state p. This 
result has been obtained by a dynamical/ variational approach in 2 , a rigorous 
proof is given in The precise statement is the following. 

Theorem 3.2. For each profile 7 : [0, 1] — > [0, 1], 

limsuplimsup^logi^^^{7r^ G Veij)} < , 
hm inf hm inf 1 log <^ {tt^ G V, (7) } > -.^(7) , 



wh' 



ere 



m . log 2i^> + 11 - 7(«)1 log ^ + (3,2) 



DYNAMICAL LARGE DEVIATIONS IN STOCHASTIC LATTICE GASES 



7 



and F G C^([0,1]) is the unique increasing solution of the non linear boundary 
value problem 

^"-(^-^)F(r^' (3.3) 
F(0) = a , F{1) = /3 . 

It is interesting to compare the large deviation properties of the stationary state 
fi with the one of /i^^, the product measure on Ejy which has the same marginals 
as i^^^, i.e. 

I^Z^IV ■■ ri{x) = 1} = p^ix) , 
where is given by (|2.2|l . It is not difficult to show that in this case 

limsuplimsup^log/i^fl{7r^ £ K(7)} < -^oil) , 

1 • .N l_N 



liminflmiinf-logAi^.^{^^eFe(7)} > -^0(7), 

£^0 Af— too 1\ 



where 



^0(7)= / I7N log ^ + [l-7(")] logy 
Jo ^ PW 1 



p{u) 1 - p{u) J 

and p is given in (|3.1(l . Notice that the functional JFq is local while J- is not. 
Moreover, it is not difficult to show [3 E| that J^) < J'- Therefore, fluctuations 
have less probability for the stationary state v^p than for the product measure 
Pa 0- This bound reflects at the large deviations level the negative correlations 
observed in ^Z.'A\ . 



4. Hydrodynamics and dynamical large deviations of the density 

We discuss the asymptotic behavior, as — > 00, of the evolution of the empirical 
density. Denote by {r}^ : i > 0} the Markov process with generator Lm and 
— TT^ {rjf). Fix a profile 7 : [0, 1] [0, 1] and assume that Tr,^^ converges to 
^{u)du as N 00. Observing the time evolution of the process, we expect 7rf to 
relax to the stationary profile p{u)du according to some trajectory pt{u)du. This 
result is usually referred to as the hydrodynamic limit. For the boundary driven 
simple exclusion process it is stated in Theorem 14 . II below jl6[ll7| . 

Fix T > and denote, respectively, by D{[Q,T],M+), £)([0, T], Eat) the space 
of A^_(--valued, EAr-valued cadlag functions endowed with the Skorohod topology. 
For a configuration rj^ in S^v, denote by P^w the probability on the path space 
D([0,r], Eat) induced by the initial state 77^ and the Markov dynamics associated 
to the generator Ljv. 

Theorem 4.1. Fix a profile 7 : [0, 1] [Oi 1] o^c^ o, sequence of configurations 
ri^ such that tt^ {rj^) converges to ^{u)du, as N ] 00. Then, for each t > 0, tt^ 
converges in Vj^n -probability to pt{u)du as N ^ 00. Here pt{u) is the solution of the 
heat equation 

dtptiu) = (l/2)Apt(u) , 

Po{u)=-i{u), ue(0,l) (4.1) 
pm = a, pt(l) = /3. 
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In other words, for each S^T > and each continuous function H : [0, 1] M we 
have 

lim P «( sup \{tt^,H) - {pt,H)\ >S) =0 

N^co ^t£[0,T] ' 

Equation (|4.1|1 describes the relaxation path from 7 to p since pt converges to 
the stationary path p as i t 00. To examine the fluctuations paths, we need first to 
describe the large deviations of the trajectories in a fixed time interval. This result 
requires some notation. 

Fix a profile 7 bounded away from and 1 : for some (5 > we have (5 < 7 < 1 — (5 
du-a.e. Denote by the following subset of D([0, T], A trajectory ttj, 
t G [0, r] is in if it is continuous and, for any t e [0, T], we have nti^du) = Xt(u)du 
for some density Xt{u) e [0,1] which satisfies the boundary conditions Aq = 7, 
At(0) = a, Xt{l) — (3- The latter are to be understood in the sense that, for each 

te [o,r], 

lim— / du \t{u) = a , 1™ T / du Xt{u) — (3 . 



sio S Jn sio S 



-s 



We define a functional /[o^T] ("It) on D{[0,T], by setting /[q^t] ('"'It) = +00 if 
TV ^ Cj and by a variational expression for tt G C^. Referring to 3, Eq. (2.4)-(2.5)] 
for the precise definition, here we note that if TTt{du) = Xt(u)du for some smooth 
density A we have 

^[o,T](^l7) = I [ dt [ dux{Xt{u))[VHt{u)]' . (4.2) 
^ Jo Jo 

Here, x(o) = a(l — a) is the mobility in the symmetric simple exclusion process and 
Ht is the unique solution of 

dtXtiu) - {l/2)AXt{u) - y[x{Xt{u))yHtiu)] , M e (0, 1) (4.3) 

with the boundary conditions Ht{0) = Ht{l) = for any t G [0,r]. As before A 
stands for and V stands for Hence, to compute /[o,t] (""It)! we first solve 
equation H4.3|l in H and then plug it in (|4.2|l . 

The following theorem states the dynamical large deviation principle for the 
boundary driven simple exclusion process. It is proven in by developing the 
techniques introduced in |14[ I23| . 

Theorem 4.2. Fix T > and a profile 7 bounded away from and 1. Consider a 
sequence rj^ of configurations associated to 7 in the sense that tt^ {rj^) converges 
to ^{u)du as N 1 00. Fix n in I?([0, T], A^_|_) and a neighborhood V^(7r) of tt of 
radius e. Then 

limsuplimsup^logP^«{7r^ G V;(7r)} < -/[o.T](7r|7) , 

liminf liminf — logP^N {tt^ G T4(7r)} > — /[o,t] ("""It) ■ 

We may now formulate the following exit problem. Fix a profile 7 and a path 
TT such that ttq = pdu, tt^ — ^du. The functional I^q j'\{'^\p) measures the cost of 
observing the path tt. Therefore, 

inf /[o,T](7r|p) 

■KT—'^du 
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measures the cost of joining p to 7 in the time interval [0, T] and 

y{l) - inf inf I[o,t]{Ap) (4-4) 

measures the cost of observing 7 starting from the stationary profile p. The func- 
tional V is called the quasi-potential. 

It is expected in general that the quasi-potential equals the rate functional of 
a large deviations principle for the empirical density under the stationary state 
^a,i3- We will see in Sectional that this is indeed the case for the boundary driven 
symmetric simple exclusion process. In particular, the quasi-potential V in (|4.4|) is 
equal to the free energy T in (|3.2|l . 

5. Dynamical approach to stationary large deviations 

In this section we characterize the optimal path for the variational problem 
(|4.4(l . As a byproduct we show that V , as defined in H4.4|l . equals as defined 
in (|3.2|l . Unless explicitly stated, the arguments presented in this section hold 
for interacting particle systems under general assumptions. We however present - 
informally - the relevant statements and proofs in the context of the one-dimensional 
boundary driven symmetric exclusion process. To simplify the notation, given a 
density path tt G D([0, T]; A^+) such that TTt is absolutely continuous with respect 
to the Lebesgue measure for each t G [0,T], TTt{du) = Xt{u)du, we shall write 
-^[o,T](A|7) for /[o,T](7r|7)- 

5.1. The reversible case. Let tptiu) be the optimal path for the variational prob- 
lem (|4.4() on the interval (— 00, 0] instead of [0, 00). In the reversible case, a = f3, 
from Onsager-Machlup we expect that it is equal to the time reversal of the re- 
laxation trajectory pt{u) solution of (|4.1|l . ipt{u) — p-t{u). We show that this is 
indeed the case. 

The cost of the path ip is not difficult to compute. By definition of Lp and by 
(Oil , dtft = -(l/2)A(p4. In particular, VHt = {V^t)lx{ft) solves so that 

Ii-oo,Q]WP) = o / dt du ■ 

Let R{a) — log{a/l — a}. Rewrite the integrand as VR{pt)'^'Pt = V{i?((pt) — 
i?(p)}V(/5t, because p is constant in the reversible case, and integrate by parts in 
space to obtain that 

/(-oo,o]('/'|p) ^ -I f [ du A^t{R{pt) ~ R{p)} ■ 







Since dt^p — — (l/2)A(pt, 5Tq{p) / 5p = R{p) — R{p), the previous expression can 
be rewritten as 

dt I dupt^^^f^^ = [ dt-^Ta{pf) = To{po) - To{p-oo) = J^oil) 
3 Jo c"^t J-oo at 

because J-o{p) — 0. This proves that V < Tq. 

The proof of LemmalObelow. with V{i?(At)-i?(/o)} instead of V {SW (Xt) / 5 Xt} , 
shows that the cost of any trajectory Af joining p to a profile 7 in the time interval 
[0,T] is greater or equal to .^0(7): 

^[0,T](A|p) > Toil) ■ 

In particular, the trajectory p is optimal and ^^(7) = J^oil)- 
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5.2. The Hamilton-Jacobi equation. We have seen in Subsection 15. II that the 
optimal path for reversible systems is the relaxation path reversed in time. In the 
non reversible case, the problem is much more difhcult and, in general, we do not 
expect to find the solution in a closed form. We first derive a Hamilton-Jacobi 
equation for the quasi-potential by interpreting the large deviation rate functional 
I[o,T]i'\p) as an action functional 

I[o,T]{\\p) = ^ r dt [\u—^--{v-'[dMu)~ (1/2) AXtiu)]}' 
^ Jq Jq X[-^t[u)) I- J 

T 

dtC{Xt,Xt) ■ 



The quasi-potential V may therefore be written as 



V(7) = inf inf / dtCiXuXt). (5.1) 

T>0 \o=P Jq 
\T='y 

From this variational formula, taking the Legendre transform of the Lagrangian, 
we derive the Hamilton-Jacobi equation for the quasi-potential: 

and 5V{'^)/5^ vanishes at the boundary. 

One is tempted to solve the Hamilton-Jacobi to find the quasi-potential and then 
to look for a trajectory whose cost is given by the quasi-potential. The problem 
is not that simple, however, because the theory of infinite dimensional Hamilton- 
Jacobi equations is not well established. Moreover, as well known, even in finite 
dimension the solution may develop caustics in correspondence to the Lagrangian 
singularities of the unstable manifold associated to the stationary solution p, see 
e.g. j20|- Finally, the Hamilton-Jacobi equation has more than one solution. In 
particular, even if one is able to exhibit a solution, one still needs to show that the 
candidate solves the variational problem H5.1|l . 

The next lemma shows that a solution W of the Hamilton-Jacobi equation is 
always smaller or equal than the quasi-potential: 

Lemma 5.1. Let W be a solution of the Hamilton-Jacobi equation H5.2|l . Then, 
W(7) — W{p) < 1^(7) for all profiles 7. 

Sketch of the proof. Fix T > 0, a profile 7, and consider a path A in Cp such 
that At = 7. We need to show that I\o,t]Wp) > W{"/) — W{p). 
The functional /[o.T](A|p) can be rewritten as 

Since SW{Xt)/SXt vanishes at the boundary, an integration by parts gives that the 
second integral is equal to 

i-T 
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Since is a solution of the Hamilton- Jacobi equation, the third integral is equal 
to 

Summing this two expressions and keeping in mind that Ht solves H4.3|l . we obtain 
that I[q j']{Mp) is greater than or equal to 

'''^di(^^^,At) = W{\t)-W{\o) = W{-i)-W[p). 

This proves the lemma. □ 
To get an identity in the previous lemma, we need the first term in H5.3|l to 
vanish. This corresponds to have VHt = V5V{\t)l5\t, i.e. to find a path A which 
is the solution of 



9tAt = (l/2)AAt-v{x(At)V- 



5\t J 

Its time reversal ipt = ^~t, t £ [— T, 0] solves 



dt^t = -(l/2)AV^t + v{x(^0 



V- 



(5.4) 



S^Pt J ' 

ip-T = 7 , 

Vt(0) = a, ^t(l)=/3. 

As we argue in the next subsection, equation (|5.4|) corresponds to the hydrody- 
namic limit of the empirical density under the time reversed dynamics; this is the 
Markov process on Ejv whose generator is the adjoint to in L2{Tjm , i^a p)- 

The next lemma shows that a weakly lower semi-continuous solution W of the 
Hamilton- Jacobi equation is an upper bound for the quasi-potential V if one can 
prove that the solution of (|5.4|) relax to the stationary profile p. 

Lemma 5.2. Let W be a solution of the Hamilton- Jacobi equation H5.2|l . lower 
semi- continuous for the weak topology. Fix a profile 7. Let ipt be the solution of 
(|5.4|) with V replaced by W. Ifipo converges p forT | 00, then V{'y) < W{j) — W{p). 

Sketch of the proof. To prove the lemma, given e > 0, it is enough to find > 
and a path (ft such that (po = p, ipr^ — 7, /[q^t^] (v'Ip) < ^^(7) ~ ^{p) + £• Fix 
T > and let ipt be the solution of equation (|5.4I) in the time interval [—T, —1] with 
initial condition ip-T — 7. Consider then an appropriate interpolation between ip^i 
and p which we again denote ipt, t £ [^1, 0]. Let ipt = "0-*: which is defined in the 
time interval [0, T]. By definition of I[o,t], 

I[a,T]{f\p) = I[OA]{v\p) + I[i,T]{'P\tp-i) ■ 

Since tp-i converges to p as T f 00, the first term can be made as small as we want 
by taking T large. The second one, by definition of ipt and by the computations per- 
formed in the proof of Lemma|0] is equal to W^(7) — VF(^_i). Since -0-1 converges 
to p and since W is lower semi-continuous we have W{p) < liminf 7-^00 W{ip-i). 
Hence limsupjn_^oQ I[o^t]{'p\p) < — W{p). This proves the lemma. □ 

Putting together the two previous lemmata, we get the following statement. 

Theorem 5.3. Let W be a solution of the Hamilton- Jacobi equation, lower semi- 
continuous for the weak topology. Suppose that the solution ipt of (|5.4I) . with V 
replaced by W, is such that 'ipo converges to p as T \ 00 for any initial profile 
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7. Then V{'j) — W{'j) — W{p). Moreover, ipt = ip-t is the optimal path for the 
variational problem (|5.1(l defined in the interval (—00, 0] instead of [0, 00). 

5.3. Adjoint hydrodynamic equation. We have just seen that equation (|5.4|) 
plays an important role in the derivation of the quasi-potential. We show in this 
subsection that (|5.4II describes in fact the evolution of the density profile under the 
adjoint dynamics. 

Consider a diff'usive interacting particle system rj^ satisfying the following as- 
sumptions. The limiting evolution of the empirical density is described by a differ- 
ential equation 

dtp = V{p) , 

where T) is a differential operator. In the symmetric simple exclusion process 'D(p) = 
(1/2) Ap. Denote by — rj^^ the time-reversed process. The limiting evolution of 
its empirical density is also described by a differential equation 

dtP^V*{p), (5.5) 

for some integro-differential operator 2?*. Moreover the empirical densities satisfy 
a dynamical large deviations principle with rate functions 

), i = l,2 

where "D^ = V and = V* for the original and time-reversed processes, respec- 
tively. In E] it is shown that 

V{p) + V*{p) = v(x(p)V^) . (5.6) 

In this general context, equation 15.411 takes the form 

dtp = -P(p) + v(x(p)v— ) . 

Therefore, under the above assumptions on the dynamics, the solution of (|5.4|l 
represents the hydrodynamic limit of the empirical density under the adjoint dy- 
namics. In particular, for non reversible systems, the typical path which creates a 
fluctuation is the time-reverse of the relaxation path of the macroscopic dynamics 
of This principle extends the Onsager-Machlup theory to irreversible systems. 



dt 



V-'{dtXt-V\Xt)) 



5.4. Optimal trajectory for the simple exclusion process. While all the ar- 
guments presented above are general, in this subsection we obtain a more explicit 
description of the optimal trajectory for the variational problem (|5.1|l for the bound- 
ary driven simple exclusion process. As a corollary we show that the quasi-potential 
is given by the expression (|3.2|l . 

Let us first show how, in this case, it is possible to obtain a solution of the 
Hamilton- Jacobi equation H5.2|l . We look for a solution of the form 

^log , ; X -0(m;p) (5.7) 



Sp{u) 1 ~ p{u) 

for some functional 0(u; p) to be determined, satisfying the boundary conditions 
m = log[a/(l - a)], m = log[/3/(l - /?)]. 
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Inserting H5.7|l into H5.2|l . we get, note that p — e'^/(l + e"^) vanishes at the 
boundary, 



= -(V log 



1-p 

(Vp,V0) + (p(l-p),(V0)2) 

e 



1 



1 



1 + e 

2 



1 



p(V(/.) 



We obtain a solution of the Hamilton-Jacobi equation if we solve the following 
ordinary differential equation which relates the functional = </'(m; p) to p 

A(t){u) 1 



1 



p{u) tie (0,1) 



(5.8) 



log 



V0 



^ m = log[a/(l - a)], (/.(I) - log[/3/(l - /3)] 
A computation shows that the derivative of the functional 

W{p)^ f du|plogp+(l-p)log(l-p) + (l-p)(/)-log(l + e' , , 

(5.9) 

is given by H5.7|l when (j){u; p) solves (I5.8|l . We note that by the change of variable 
4> = log[F/(l - F)] equation (|5.8|) becomes H3.3|l and (|5.9|l becomes H3.2|l . 

We next show that satisfies the hypotheses of Theorem 15 .31 In this case (|5.4f) 
(with V replaced by W), after shifting the time interval from [— T, 0] to [0, T], is 
given by the equation non local in space 



-Ap-V{p[l-p]V(/)(p)} tie (0,1) 



(5.10) 



/9t(0) = a, Pt(l)=/3 
. Po{u) = 7(u) 

where (/)(jo) is to be obtained from p by solving H5.8(l . Since 0(p) = log[p/(l — p)], 
we see that p is also a stationary solution of l|5.10|) . 

Equation (|5.10() can be related to the heat equation as follows. Let pt be the 
solution of H5.10|l and introduce F = Ft{u) as 

g0(";pt) 

FAu) = -7 . (5.11) 

it is not too difficult, see [3 Appendix B], to check that Ft{u) satisfies the heat 
equation 

1 



dtFtiu) ^ -/^Ft{u) UG(0,1) 
Ft{Q) = a, Ft{l) = (3, 

Fo{u) 



(5.12) 



1 + e0(";7) 

Conversely, given Ft{u) which solves H5.12|l . by setting 



pt{u) ^ Ftiu) + Ft{u)[l - Ft{u)] 



AFtju) 
[VF,{u)Y 
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a computation shows that pt solves (jS.lOII . 

We have thus shown how a solution of the (non local, non linear) equation 
l)5.1U|l can be obtained from the linear heat equation by performing the non local 
transformation H5.11|l on the initial datum. In particular, since the solution Ft{u) 
of (|5.12() converges as t ^ cxd to p, we see that the functional W{p) given in (|5.9() 
satisfies the hypotheses of Theorem 15.31 Since W{p) = we thus have that the 
quasi-potential of the boundary driven simple exclusion process is given by the 
functional 1)5. 9|) where ip solves (|5.1U|I . 

To conclude the dynamical proof of Theorem 13.21 we have to identify the rate 
function J-' of the invariant measure with the quasi-potential V. 

5.5. Free energy and quasi-potential. Bodincau and Giacomin j]J]j, adapting 
to this infinite dimensional setting the method introduced by Freidlin and Wentzell 
|18| in the context of small perturbations of dynamical systems, proved the following 
theorem which identifies V and T. 

Theorem 5.4. Let I^q t] be the rate function in Theorem \4.^ and define the quasi- 
potential as in (14.411 . Then the empirical density under the stationary state satisfies 
a large deviation principle with rate functional given by the quasi-potential. 

An explicit description of the quasi-potential, as the one here discussed for the 
boundary driven simple exclusion process, is not always possible. There are other 
few one dimensional boundary driven models for which a similar representation has 
been obtained. This class includes the weakly asymmetric exclusion processes, the 
zero range processes and the KMP model |P |7| . 



We examine in this section the current fluctuations over a fixed macroscopic time 
interval. In particular we discuss the law of large numbers and the dynamical large 
deviations principle for the empirical current. We state these results in the context 
of the boundary driven symmetric exclusion process but similar results hold for 
more general dynamics and for periodic boundary conditions. 

Consider the boundary driven symmetric simple exclusion process defined in 
Section 121 For < a; < - 1, denote by jx,x+i the rate at which a particle jumps 
from X io X + 1 minus the rate at which a particle jumps from a; + 1 to a;. For 
a; = 0, this is the rate at which a particle is created minus the rate at which a 
particle leaves the system. A similar interpretation holds at the right boundary. 
An elementary computation shows that 



In view of (|2.2() . under the invariant measure the average of jx,x+i is 



Given a bond {x, a:-}-!}, < a; < A^— 1, let j^^'^+i (resp. J^'^^'^) be the number of 
particles that have jumped from x to (resp. a; -1-1 to a;) in the time interval [0, t]. 
Here we adopt the convention that J^'^ is the number of particles created at 1 and 
that J^'^ represents the number of particles that left the system from 1. A similar 
convention is adopted at the right boundary. The difference W^''^'^^ = — 



6. Asymptotic behavior of the empirical current 
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is the net number of particles flown across the bond {x,x + 1} in the time 
interval [0,t]. Let us consider the stationary process Pi/™^; i-e- the boundary driven 
symmetric simple exclusion process in which the initial condition is distributed 
according to the invariant measure J^^^- A simple martingale computation shows 



N> 1, x^O 



that W^"^'^^/{Nt) converges, ast ^ oo, to (a — /3) in probability. Namely, for each 



, TV — 1, and 5 > 0, we have 



lim I 

t — '■OO 



Wi 



x,x-\-l 



Nt 



- (a - /?) 



> S 



= 



Let j\4 be the space of bounded signed measures on [0, 1] endowed with the weak 
topology. For t > 0, define the empirical integrated current W^^ 6 as the finite 
signed measure on [0, 1] induced by the net flow of particles in the time interval 
[0,t]: 

N-l 



N- 



^x/N 



x=0 



Notice the extra factor iV^^ in the normalizing constant which corresponds to the 
diffusive rescaling of time. In particular, for a function F in C([0, 1]), the integral 
of F with respect to , also denoted by {W^^ , F), is given by 



N-l 

x=0 



F{x/N)Wt 



X,X~\-1 



(6.1) 



It is not difficult to prove the law of large numbers for the empirical current starting 
from an initial configuration associated to a density profile. 

Proposition 6.1. Fix a profile 7 and consider a sequence of configurations rj^ 
such that TT^ (rj^) converges to j{u)du, as N 1 00. Let p be the solution of the heat 
equation H4.1|l . Then, for each T > 0, S > and F in C([0, 1]), 

i-T rl 



lim 



{W^\F) + (1/2) 



dt 



F{u)Vpt{u) du 



> 5 







This result states that the empirical current Wl converges to the time integral 
of — (l/2)Vpt(u), which is the instantaneous current associated to the profile pt. 
Thus, if we denote by ^(7) = — (1/2)V7 the instantaneous current of a profile 7, 
we have that 



lim = [ ds 



wiPs) 



in probability. Proposition lft.ll is easy to understand. The local conservation of the 
number of particles is expressed by 



Vt{x) - r]a{x) = Wt 



w; 



x,x+l 



It gives the following continuity equation for the empirical density and current. 
Let G : [0, 1] ^ M be a smooth function vanishing at the boundary and let 

{VnG){x/N) = N{G{x + 1/N) - G{x/N)}. Then, 

(7rf,G)-(^o^,G) = (W^f,V^G). 

The previous identity shows that the empirical density at time t can be recovered 
from the initial state and the empirical current at time t. In contrast, the empirical 
density at time t and at time determines the empirical current at time t only up 
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to a constant. Letting "f oo in the previous identity, since tt converges to the 
sohition of the heat equation 14.1|l . an integration by parts gives that 

{Wt,VG) = (pt,G)-(po,G) - ^j\s{Ap,,G) = J^s {Vp,,^G) . 

where Wt is the hmit of . 

After proving this law of large numbers for the current, we examine its large 
deviations properties. To state a large deviations principle for the current we need to 
introduce some notation. Fix T > and recall that we denote by 10(7) = — (1/2)V7 
the instantaneous current associated 7. For a density profile 7 and a path W in 
D{[0,T], A4), denote by wt — Wt and let p]'^ the weak solution of 

dtpt + ywt = , 

Po{u) = 7(u) , (6.2) 
Pt(0) = a, pt{l)^f3. 

We note that the trajectory p]'^ is the one followed by the density profile if the ini- 
tial condition is 7 and the instantaneous current is w. As for the empirical density, 
the rate functional for the empirical current is given by a variational expression. 
Referring to W for the precise definition, we here note that for trajectories W in 
£'([0, T],Ai) the rate functional is finite only if the associated density path p]'' du 
belongs to C([0, T], A^+); moreover when Vl^ is a smooth path we have 

%T](W^|7) = \ f dt{—^ {Wt~w{pT'^)Y) . (6.3) 

The following theorem is proven in [Hj in the case of periodic boundary condition. 
The proof can easily be modified to cover the present setting of the boundary driven 
simple exclusion process. 

Theorem 6.2. Fix T > and a smooth profile 7 bounded away from and 1. 
Consider a sequence ry^ of configurations associated to 7 in the sense that tt^ {rj^) 
converges to j{u)du as N '] 00. Fix W in D{[0,T], Jvl) and an associated neighbor- 
hood V!^{W) of radius e. Then, 

limsuplimsup4logP^iv{W^^ e ^^(W^)} < -^[o,t]{W\j) , 

liminf Irniinf ^logP,,iv{T4^^ e ye(VK)} > -Z[o,t] (VKI7) . 

Since the trajectory of the empirical density can be recovered from the evolu- 
tion of the current and the initial condition, the large deviations principle for the 
empirical density stated in Theorem 14 . 21 can be obtained from the large deviations 
principle for the current by the contraction principle, see jBj for the proof. 

7. Large deviations of the time averaged empirical current 

In this section we investigate the large deviations properties of the mean empiri- 
cal current W^/T as we let first — > 00 and then T — > 00. As before, unless stated 
explicitly, the analysis carried out in this section does not depend on the details of 
the symmetric simple exclusion process so that it holds in a general setting. 

Since the density is bounded, for T large the time averaged empirical current 
must be constant with respect to the space variable u. This holds in the present 
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one-dimensional setting; in higher dimensions the condition required would be the 
vanishing of the divergence. Indeed, if this condition were not satisfied we would 
have an unbounded (either positive or negative) accumulation of particles. We next 
discuss the asymptotic probability that the time averaged empirical current equals 
some fixed constant. 

For a smooth profile 7 bounded away from and 1, let $ : — > [0, +00] be the 
functional defined by 

~ f inf — inf Tin t](W\j) if J(du) — qdu for some o G M 

1+00 otherwise 
where AT,q stands for the set of currents with time average equal to q 

AT,q.= [w ^D{[0,T]-M) : - J dtWt{du)=qduy 

It is not difficult to show that $ is convex. In the present context of the boundary 
driven simple exclusion process, it is also easy to verify that the functional $ does 
not depend on on the initial condition 7. We emphasize however that, in the 
case of periodic boundary condition $ depends on 7 only through its total mass 
Jduj{u). Indeed, we may start by driving the empirical density from a profile 7 to 
a profile 7' in the time interval [0, 1] paying a finite price, note that in the periodic 
case 7 and 7' must have the same mass. As T f c», this initial cost vanishes 
and the problem is reduced to the original one starting from the profile 7'. Let 
us finally introduce $ as the lower semi-continuous envelope of i.e. the largest 
lower semi-continuous function below $. The next theorem states that, as we let 
first iV t cx) and then T f cxo the time averaged empirical current /T satisfies a 
large deviation principle with rate function We refer to |S] for the proof which 
is carried out by analyzing the variational problem infive^T,,, ^~^^[o,T](W^l7) ^ 
T I 00 and showing that it converges, in a suitable sense, to the variational problem 
defining $. 

Theorem 7.1. Fix T > and a smooth profile 7 hounded away from and 1. 
Consider a sequence i]^ of configurations associated to 7 in the sense that tt^ (rj^) 
converges to ^{u)du as N 00. Fix J ^ A4 and a neighborhood V^(J) of radius e. 
Then, 



lim sup lim sup lim sup — — log I 



lim inf lim inf lim inf — — log 

£^0 T^oo Af^co TN 



1 . „ ri 



^W^ e V,{J) > -$(J) 



A result analogous to Theorem 17. II can be proven for other diffusive interacting 
particle systems. Consider a system with a weak external field E = E{u), whose 
hydrodynamic equation, describing the evolution of the empirical density on the 
macroscopic scale, has the form 

dtpt = V{D{pt)Vpt) - V{x{pt)E) . (7.2) 

where D{p) is the diffusion coefficient and xip) is the mobility. For the symmetric 
simple exclusion process D = 1/2 and x{p) = p(l ~ p)- In the general case, 
the large deviations functional T[o,t]('|7) has the same form (|6.3ll with ^(7) = 
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— I?(7)V7 + x{l)E ^-nd p^'^ the solution of H6.2|l . For systems with periodic 
boundary conditions, the boundary conditions in Il6.2(l is modified accordingly. In 
the remaining part of this section we analyze the variational problem 1)7. l|l for 
different systems and show that different scenarios are possible. 

A possible strategy for minimizing ^[q^t] {"u^ll) with the constraint that w e •^T,q, 
i.e. that the time average of w is fixed, consists in driving the empirical density to a 
density profile 7*, remaining there most the time and forcing the associated current 
to be equal to q. This is the strategy originally proposed by Bodineau and Derrida 
12]. In view of H().3|) the asymptotic cost, as T f 00, of this strategy is 

i( [g + W)V7*] , ^ [<Z + i?(7*)V7* 

If we minimize this quantity over all profiles 7* we obtain a functional U which 
gives the cost of keeping a current q at a fixed density profile: 

U(q) := inf i( [q + Dip)\7p] , ^ [-7 + D{p)Vp] ) . (7.3) 

where the infimum is carried out over all smooth density profiles p = p{u) bounded 
away from and 1 which satisfy the boundary conditions p{0) = a, p{l) — j3. 
As observed above, for boundary driven systems all density profiles are allowed 
while for periodic boundary condition only profiles with the same total mass m = 
J^dup{u) are allowed. In the latter case, the functional U depends on the total 
mass m and is denoted by [/,„. 

As proven in 0] El |2| , for the symmetric simple exclusion process the strategy 
above is the optimal one, i.e. ^ — U. It is in fact not difficult to show that in 
this case U is lower semi-continuous, so that <& = $. More generally we have the 
following result. 

Lemma 7.2. Let E = 0. IJ D{p)x"{p) < D'{p)x'{p) for any p, then $ = JJ. 

Besides the symmetric simple exclusion process, the hypothesis of the lemma 
is also satisfied for the zero range model, where D{p) — '^'{p) and x{p) ~ ^(p) 
for some strictly increasing function 4" : K_|_ — > K_|_, and for the non interacting 
Ginzburg-Landau model, where p S R, D(p) is an arbitrary strictly positive func- 
tion and x{p) is constant. 

For systems with periodic boundary condition we have shown |S] that the 
profile which minimizes Um is the constant profile if l/x{p) is convex. 

Lemma 7.3. Let E — 0. If the function p > 1/xip) is convex, then 

1 fl2 



2x(m) 

and the constant profile p{u) — ni is optimal for the variational problem (|7.3|l . 

The assumption of this lemma is satisfied by the symmetric simple exclusion 
process as well as by the KMP model |3|22j|, where D{p) = 1 and x{p) = ■ 

As first discussed in jj, the above strategy is not always the optimal one, i.e. 
there are systems for which $ < [/. In we interpreted this strict inequality as 
a dynamical phase transition. In such a case the minimizers for 1)7. l|l become in fact 
time dependent and the invariance under time shifts is broken. We now illustrate 
how different behaviors of the variational problem 1)7.1)1 leads to different dynamical 
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regimes. We consider the system in the ensemble defined by conditioning on the 
event {T)~^Wr^ {du) — qdu, q E M., with N and T large. The parameter q plays 
therefore the role of an intensive thermodynamic variable and the convexity of $ 
expresses a stability property with respect to variations of q. 

If ^{q) — U{q) and the minimum for (|7.3|l is attained for p — p{q) we have a state 
analogous to a unique phase: by observing the system at any fixed time t = 0{T) 
we see, with probability converging to one as N,T ^ oo, the density tt^ ^ p{q) 
and the instantaneous current Wf^ ^ q. 

While the functional $ is always convex, U may be not convex; an example of 
a system with this property is given in 0. If $ is equal to the convex envelope 
of U, we have a state analogous to a phase coexistence. Suppose for example 
<? = Ml + (1 -P)92 and U{q) > U**{q) = pU{qi) + {l-p)U{q2) for some p, 171,92; 
here U** denotes the convex envelope of U . The values p, qi, 52 are determined by 
q and U. The density profile is then not determined, but rather we observe with 
probability p the profile p{qi) and with probability 1 — p the profile p{q2)- 

Consider now the case in which a minimizer for H7.1(l is a current path wt not 
constant in t and suppose that it is periodic with period r. We denote by pt the 
corresponding density. Of course we have Jg dtwt = q. In this situation we 
have in fact a one parameter family of minimizers which are obtained by a time 
shift a G [0,t]. This behavior is analogous to a non translation invariant state in 
equilibrium statistical mechanics, like a crystal. Finally, if the optimal path for 
(|7.1(l is time dependent and not periodic the corresponding state is analogous to a 
quasi-crystal. 

The explicit formula for Um derived in Lemma |7. 31 permits to show that under 
additional conditions on the transport coefficient D and a dynamical phase 
transition occurs. We discuss only the case of periodic boundary conditions. In 
this situation a time-averaged current q may be produced using a traveling wave 
density profile, pt(u) = poiu — vt), with velocity w ~ g. Assume now that E = 
and the function p 1-^ x{p) is strictly convex for p = m. Then, for sufficiently large 
q, the traveling wave strategy is more convenient than the one using the constant 
profile TO 01 El . In particular, if p i— > l/)({p) is convex so that Lemma |7. 31 can be 
applied, we have 

< Um{q) (7.4) 

for sufficiently large q. In the KMP model the above hypotheses are satisfied for 
any to > 0; we can thus conclude that a dynamical phase transition takes place for 
sufficiently large time-averaged currents. 

The above analysis can also be applied to the weakly asymmetric simple exclusion 
process 0. It yields that if \E/q\ > [to(1 — m)]~^ for q large there exists a traveling 
wave whose cost is strictly less than the one of the constant profile p{u) = to. 
The analysis in 9 suggest however that the strict inequality 1)7. 4|l holds also in 
this case. Moreover, the numerical simulations in |2] indicate the existence, for the 
weakly asymmetric simple exclusion process, of a critical current q* below which 
the optimal profile is constant and above which the optimal profile is a traveling 
wave. 
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